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ABSTRACT 


Laser-induced  ablation  of  fibre-resin  composites  is  considered  using 
a  steady-state,  one-dimensional  model.  In  the  analysis,  it  is  attempted  to 
clarify  the  physical  and  chemical  aspects  of  ablation  rather  than  to  predict 
ablation  rates.  It  is  shown  how  suitable  measurements  may  be  used  in 
conjunction  with  this  model  to  determine  loss  rates  for  combined  gaseous  and 
combined  non-gaseous  products  of  ablation,  and  to  determine  the  average 
enthalpy  of  each  of  these  two  losses.  The  model  is  restricted  to  incident 
irradiances  *ti  excess  of  approximately  30  MW  and  to  cases  in  which  the 
ablation  rate  varies  with  irradiance  in  a  linear  fashion.  Because  some  of 
the  required  measurements  may  be  difficult  to  make  accurately,  this  analysis 
should  be  regarded  as  a  preliminary  to  more  detailed  work.^ 
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described  here*  This  model  is  restricted  to  linear  variations  of  ablation 
rate  with  incident  laser  irradiances,  and  to  irradiances  above  approximately 
30  MW  m~^,  for  which  some  simplifying  assumptions,  not  otherwise  valid,  can  be 
made.  In  addition,  only  steady-state  ablation  is  considered.  Shortcomings 
in  the  model  and  in  its  application  are  discussed* 


2.  PHYSICAL  MODEL 


Figure  1  shows  the  model  adopted  here  for  the  steady-state  ablation 
of  a  fibre-resin  composite*  Physically,  it  is  based  on  standard  concepts 
[1-4].  The  model  is  one-dimensional  in  the  x  direction  with  x  =  0  fixed  on 
the  ablating  surface  so  that  material  appears  to  move  in  the  positive  x 
direction  and  alters  composition  chemically  as  it  moves*  Three  zones  of 
composition  exist  which  have  fixed  boundaries  for  a  given  ablation  rate  v. 


In  zone  A,  virgin  material  extends  effectively  to  x  =  -*•  At  an 
unspecified  value  of  x,  zone  B  begins  in  which  some  decomposition  of  the  resin 
occurs  at  a  rate  that  increases  with  the  temperature  T(x)  according  to  the 
Arrhenius  law  [8],  However,  the  resin  pyrolysis  is  not  advanced  to  the  stage 
in  which  a  hard  char  is  formed*  The  model  assumes  that,  in  this  zone,  no 
decomposition  of  the  fibres  occurs.  Thus  fibres  like  kevlar,  which  have 
similar  degradation  temperatures  to  resins,  are  not  included  in  the  model. 

Zone  C  is  the  surface  char  zone,  and  is  assumed  to  begin  at  a  measurable 
x  =  x^ .  Resin  pyrolysis  in  this  zone  produces  a  hard  char,  and  decomposition 
or  melting  of  the  fibre  also  occurs.  The  fact  that  surface  char  layers 
usually  have  a  reasonably  well-defined  thickness  is  well  established  [1-3], 
and  is  a  consequence  of  the  large  difference  in  the  temperature  stability  of 
char  compared  with  the  stability  of  undegraded  resins.  In  zone  C,  gas 
produced  from  the  resin  in  zone  B  may  undergo  further  decomposition  before 
emerging  at  x  =  0,  Much  of  the  surface  char  will  be  lost  as  gas,  but  some  of 
it,  molten  material,  and  undecomposed  fibres  may  be  thermomechanical ly  removed 
from  positions  close  to  x  =  0.  The  analysis  here  assumes  that  the 
thermomechanical  removal  rate  cannot  be  measured  directly  by  any  means.  Any 
reaction  with  environmental  gases  occurs  only  in  the  gas  phase  external  to  the 
ablating  body  because  the  evolving  gases  convectively  block  environmental 
gases. 


A  fraction  6  of  the  incident  irradiance  PL  penetrates  the  gases 
above  the  surface.  All  of  the  5Pl  is  assumed  to  be  absorbed  at  x  =  0  because 
chars  have  very  high  absorption  coefficients  [9]  at  the  laser  wavelengths  of 
interest. 


The  profiles  of  temperature  T(x)  (above  ambient  temperature)  in 
zones  B  and  C  may  differ  but  are  assumed  to  be  linear,  with  Tfx^  =  T1  small 
compared  with  T(0)  =  Tg  at  the  surface.  (In  fact,  T(x)  in  zone  C  is  unlikely 
to  be  linear,  but  the  validity  of  this  assumption  is  defended  in  Sec.  3.3). 

In  zone  A,  T(x)  is  of  no  importance  even  for  conditions  prior  to  steady-state 
because  the  line  of  demarcation  between  zones  A  and  B  may  be  chosen  at  any  x 
value  for  which  T  <<  T1 ;  that  is,  the  enthalpy  of  the  virgin  material  may  be 
made  arbitrarily  small.  Thus  the  zone  of  virgin  material  is  included  in  the 

model  only  for  the  purpose  of  completing  the  physical  description. 


A  SEMI-ANALYTICAL  MODEL  FOR  LASER- INDUCED 


ABLATION  OF  FIBRE-RESIN  COMPOSITES 


1  *  INTRODUCTION 


For  military  and  industrial  applications,  high  energy  lasers  may  be 
used  to  penetrate  or  cut  fibre-resin  composite  materials.  In  the  case  of 
continuous-wave  lasers,  incident  irradiances  are  typically  between  30  MW  m” 
and  1  GW  m~2.  Ablation  of  fibre-resin  heat  shields  has  been  extensively 
studied  for  aerodynamic  heating  rates  below  25  MW  m”2  [ 1 —4]  and  much  of  this 
work  is  relevant  to  laser-induced  ablation,  which  has  received  less  attention 
[5-7] • 


Except  in  some  unusual  circumstances,  the  complexity  of  the  physical 
and  chemical  ablation  processes  for  fibre-resin  composites  is  such  that  no 
method  currently  exists  which  enables  the  ablation  rate  of  these  materials  to 
be  accurately  predicted.  For  laser-induced  ablation,  three  basic 
inadequacies  in  present  theoretical  analyses  have  been  identified  [5], 

Firstly,  the  loss  rate  and  enthalpy  of  solid  or  liquid  material  removed 
thermomechanically  cannot  be  calculated,  and  appears  to  be  very  difficult  to 
measure  directly.  Such  removal  may  be  caused  by  spallation  arising  from 
pressure  of  pyrolysis  gases  or  by  rapid  release  of  thermal  stresses  in  the 
char,  and  if  an  imposed  gas  flow  exists,  aerodynamic  shear  and  pressure  forces 
may  assist  this  removal.  Secondly,  thermophysical  constants  (e.g.  char 
thermal  conductivity)  and  reaction  rate  constants  are  often  poorly  known. 

These  constants  determine  in  part  the  final  gaseous  ablation  products  and 
their  relative  concentrations.  Thirdly,  the  laser  radiation  may  be  partially 
absorbed  before  reaching  the  material  by  pyrolysis  gases  or  by  products  of 
reactions  between  these  gases  and  environmental  gases  (usually  air),  Hie 
magnitude  of  this  ”gas  blockage”  effect  is  difficult  to  estimate. 

This  report  presents  a  simple  analysis  which,  in  conjunction  with 
experiments  on  a  given  material,  enables  quantitative  information  to  be 
obtained  on  some  aspects  of  the  three  factors  mentioned  above,  and  which 
attempts  to  clarify  the  physics  of  the  ablation  processes.  An  ability  to 
predict  ablation  rates  awaits  more  detailed  analyses  applied  individually  to 
each  particular  composite,  but  these  analyses  may  be  based  on  the  model 
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Some  assumptions  are  made  concerning  the  two  types  of  thermophysical 
parameters  of  relevance,  namely  mass  density  and  thermal  conductivity*  The 
resin  density  in  zone  B  and  the  density  pc  in  zone  C  are  assumed  to  vary 
with  x  at  rates  that  are  not  significantly  greater  than  the  corresponding 
profiles  for  T(x).  Many  complex  factors  influence  these  density  profiles, 
but  the  validity  of  this  assumption  is  supported  by  considerable  experimental 
evidence  [1,2]*  A  non-critical  assumption  is  that  the  thermal  conductivity 
kfi  in  zone  B  is  known  or  is  easily  measured*  This  is  a  good  assumption 
because  the  conductivity  of  the  undecomposed  fibres  is  the  dominant 
contribution  to  kfi  [10],  and  this  conductivity  is  known  or  easily  determined. 
An  important  feature  of  the  present  model  is  that  the  conductivity  in  zone  C 
is  not  assumed  to  be  known.  Several  other  analyses  assume  that  char 
conductivity  is  a  measurable  quantity,  but  this  appears  to  be  a  very  poor 
assumption  in  view  of  the  dominance  of  radiative  transport  through  porous  char 
at  high  temperatures  [11]. 


3.  ANALYSIS 


3  •  1  Theoret  ical  Approach 

The  model  utilizes  conservation  of  energy  and  continuity  of  mass 
with  the  main  aim  being  to  determine  from  measurements  the  rate  at  which 
material  is  removed  thermomechanically,  and  the  average  enthalpy  of  this 
material  at  the  time  of  removal.  Individual  knowledge  of  these  two 
quantities  provides  a  better  appreciation  of  the  ablation  process  than  simply 
a  knowledge  of  their  product.  The  ultimate  aim  of  the  model  is  to  acquire  a 
capacity  to  predict  the  behaviour  of  composites  under  various  conditions  of 
irradiance. 


Consider  first  the  hypothetical  case  in  which  all  material  is 
removed  in  gaseous  form.  Conservation  of  energy  in  the  steady-state  gives 
[12] 


6PL-PR=P-VHG  (1> 

where  PR  is  the  energy  lost  by  radiation  and  convective  cooling  per  unit  area 
per  unit  time,  is  the  density  of  the  virgin  material,  v  is  the  ablation 
rate,  and  Hq  is  the  enthalpy  of  the  gases  at  the  time  of  removal,  measured 
with  respect  to  the  virgin  material.  This  enthalpy  represents  an  upper  limit 
to  the  average  enthalpy  of  ablation  products  in  the  more  general  case  in  which 
a  fraction  of  the  material  is  removed  thermomechanically.  In  this  latter 
case,  Eq.  (1)  becomes 


6Pl  -  PR  *  p.  v[(1  -  F)HG  +  FHS]  (2) 


where  F  is  the  mass  fraction  removed  thermomechanically,  and  Hg  is  the  average 
enthalpy  of  the  non-gaseous  products  at  the  time  of  removal. 
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the  aim  here  is  to  determine  Hg  from  Eq.  (2)  after  evaluating  the 
other  parameters  from  experiments  and  from  theory  involving  known  chemical 
details.  Most  of  the  present  analysis  is  concerned  with  the  evaluation  of  F 
(Section  3.3).  Progress  towards  this  aim  is  possible  if  the  analysis  is  now 
subjected  to  two  main  restrictions. 

Firstly,  the  parameters  5,  PR,  HQ,  Hg,  and  F  are  all  considered  to 
be  constants  over  the  PL  regime  of  interest.  From  Eq.  (2)  this  means  that  v 

is  linearly  dependent  on  P^.  As  explained  in  Sec.  3.2  this  restriction  also 
means  that,  for  gaseous  products,  only  equilibrium  decomposition  is 
considered.  Uie  variation  v(PL)  can  be  determined  experimentally  [7,12],  and 
if  found  to  be  linear,  this  restriction  is  likely  to  apply.  Little 
possibility  exists  that  compensating  variations  in  HQ  and  Hg  with  PL  will 
result  in  a  linear  v(PL)  because  Hg  will  tend  to  be  constant; 
thermomechanical  removal  occurs  when  the  material  has  attained  a  particular 
physical  condition.  (If  F  is  non-constant,  this  is  observable  independently 
of  v(PL)  -  see  Sec.  3.3.)  Linear  variations  for  v(P^)  have  been  observed  for 
several  ablating  composites  [6,7].  However,  a  non-linear  v(PL)  has  been 
observed  for  kevlar-epoxy  ablating  in  an  airflow,  and  this  was  attributed  to  a 
strong  dependence  of  HG  on  PL  arising  from  nonequilibrium  decomposition 
effects  [7]. 


The  second  major  restriction  here  is  that  only  the  regime  of 
irradiance  PL  in  excess  of  approximately  30  MW  m“2  is  considered.  As 
explained  later  in  Secs.  3.2  and  3.5,  this  effectively  enables  HG  and  PR  to  be 
evaluated.  This  is  not  a  serious  restriction  because  in  practice  lasers  are 
rarely  used  at  irradiances  below  30  MW  nT2  for  cutting  or  penetrating  non- 
metals. 


3.2  Determination  of  HG 

_2 

For  PL  >  30  MW  m  ,  the  gaseous  ablation  products  may  be  assumed  to 
be  a  small  number  of  elements  or  simple  compounds  such  as  C,  Si,  CO,  SiO,  H, 
H2,  CH,  C2H  and  C^H  [5].  Relative  concentrations  of  these  product  species 
may  be  predicted  with  sufficient  accuracy  by  analytical  or  numerical  methods 
in  order  to  identify  a  small  number  that  are  dominant.  This  is  so  because 
uncertainties  in  chemical  kinetic  rate  constants  become  less  significant  at 
temperatures  sufficiently  high  to  induce  almost  complete  dissociation  of 
hydrocarbons.  Exceptions  occur,  for  example,  when  the  heating  rates 
experienced  by  a  material  during  ablation  are  faster  than  those  rates  that 
allow  equilibrium  decomposition  conditions.  Such  nonequilibrium  effects 
cause  the  gas  enthalpy  HG  to  vary  with  PL,  and  thus  v(PL)  is  non-linear  [7]. 
The  present  model  accommodates  only  those  nonequilibrium  events  which  lead  to 
non-gaseous  ablation  products,  i.e.  thermomechanically-removed  products. 


If  attention  is  restricted  to  equilibrium  decomposition  for  gaseous 
products  (i.e.  high  P^  and  linear  v(PL)),  HQ  may  be  obtained  in  part  from  the 
heat  of  formation  of  the  material  considered  to  have  been  formed  from  the  few 
dominant  gaseous  ablation  products.  Since  these  gases  are  simple  compounds 
or  elements,  their  own  heats  of  formation  are  low  and  do  not  strongly 
influence  the  overall  heat  of  formation.  Thus  considerable  uncertainty  in 
the  gas  composition  may  be  tolerated.  In  order  to  obtain  Hq,  the  sensible 


J 
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heat  content  CpTg  of  the  gas  must  be  added  to  the  overall  heat  of  formation 
(where  Cp  is  the  average  specific  heat  of  the  gases). 

3.3  Determination  of  F 


This  is  the  only  parameter  which  requires  for  its  evaluation  a 
consideration  of  the  zones  in  Fig.  1.  For  zone  C,  the  mass  flowing  in  at 
x  *  x-|  equals  that  flowing  out  at  x  =  0;  i.e. 


vpF  +  vpR1  +  Mgr  =  fMGR  +  M qq  +  Mg 


(3) 


where  pp  is  the  mass  density  of  undecomposed  fibres,  pR-j  is  the  mass  density 
of  non-gaseous  resin  remaining  at  x^ ,  MGR  is  the  mass  of  gas  evolved  through 
x1  from  the  resin  per  unit  area  per  unit  time,  f  is  the  mass  fraction  of  this 
gas  that  emerges  at  x  =  0'  mgc  is  the  mass  of  gas  evolved  through  x  *  0  from 
char  decomposition  per  unit  area  per  unit  time,  and  Mg  is  the  rate  at  which 
mass  is  lost  thermomechanically  per  unit  area.  For  zone  B,  continuity  of 
mass  gives 


VPoo  =  Vpp  +  VpR1  +  MGR 


(4) 


Equations  (3)  and  (4)  give 


fM _ 

+ 

+  M_ 

GR 

GC 

S 

fMGR 

+  M„ 
GC 

+  Fvp 

00 

(5) 

It  is  now  assumed  that  MGR  can  be  determined  from  a  simple 
measurement  using  a  sample  of  the  resin.  When  a  resin  is  heated,  a  mass 
fraction  of  it  forms  a  relatively  stable  char  while  the  remaining  fraction  r 
is  evolved  as  gas.  Measurements  of  r  are  not  uncommon,  and  a  typical  value 
for  epoxy  is  0.8  [13).  The  heating  rates  employed  in  these  measurements  are 
slower  than  those  induced  by  a  laser  because  observation  of  resin  weight  loss 
with  temperature  is  necessary  in  the  former  case.  Since  the  pyrolysis 
behaviour  at  low  heating  rates  may  not  be  the  same  as  those  at  high  heating 
rates,  it  is  assumed  here  that  the  fraction  r  is  the  same  in  both  cases,  but 
that  the  temperature  for  stable  char  formation  (T-j  in  the  case  of  laser 
heating)  may  be  different.  (Determination  of  is  discussed  later  in  this 
section).  Thus  MQR  in  Eq.  (5)  can  be  replaced  by 


mGR  3  rvPR  (6) 

where  pR  is  the  known  density  of  virgin  resin. 
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The  additional  assumption  is  now  made  that  f  *  1,  i.e*  that  although 
further  dissociation  of  the  gas  represented  by  may  occur  in  zone  C,  none 
of  this  gas  is  deposited  there.  Reactions  that  may  occur  between  this  gas 
and  char  in  zone  C  are  viewed  as  exchange  reactions  (mainly  of  carbon)  without 
loss  to  M^r.  This  non-interactive  nature  of  MQR  and  zone  C  may  be  expected 
because  for  the  very  high  temperatures  in  zone  C  (about  3,500  K  at  the 
surface),  decomposition  processes  are  dominant  over  recombination  i. rocesses 
[4,5]. 


In  order  to  determine  for  use  in  Eq.  (5),  the  Arrhenius  law  is 
assumed  for  the  decomposition  of  the  char  [8]: 

MGC  =  Xl  PC  Bc  exp(-Ec/(RT))  dx  (7) 

where  pc  *  pc<x)  *s  the  density  of  solid  plus  liquid  in  the  char,  Bc  is  the 
frequency  factor,  E^,  is  the  activation  energy,  and  R  is  the  gas  constant. 

Since  the  chemical  composition  of  zone  c  is  very  complex,  the  use  of  a  single 
expression  of  this  sort  represents  a  gross  simplification.  However,  use  of 
Eq.  (7)  is  somewhat  justified  if  it  is  viewed  not  in  terms  of  describing  a  set 
of  chemical  processes  but  as  an  expression  for  which  an  empirical  fit  is  to  be 
obtained  based  mainly  on  physical  considerations. 


An  approximate  solution  for  the  integral  term  in  Eq.  (7)  is  possible 
because  invariably  EC/(RT)  >>  1 .  Thus  the  exponential  term  increases  very 
rapidly  with  T,  while  p  varies  slowly;  i.e.  most  of  the  contribution  to  MGC 
comes  from  the  part  of  the  char  near  the  hot  surface.  Therefore,  p  (x)  may 
be  approximated  by  the  density  p  at  the  surface,  and  T  may  be  expanded  near 
Tsas 


where 


T  =  Tg  +  bx 


(8a) 


b  =  (T,  -  Ts)/Xl 


(8b) 


thus 


T-1  »  (Ts  -  bx)/[(Tg  -  bx) (Tg  +  bx)] 


=  (Tg  -  bx)/(T|  -  b2  x2) 


“  Tg1  -  bx/Tg  (since  b2  x2  <<  t|) 


(9) 
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Hence 


x  »  (Tg/b) (Tg1  -  T"1)  (to) 

and 

dx  »  -(Ig/bJdCr”1 )  (11) 


Using  Bq.(11)  to  integrate  Bq.(7)  and  putting  exp(-Ec/(RTg) )  »  exp(-Ec/(RT. ) ) 
gives 


M 


GC 


psBcRTs 

bE„ 


exp(-Ec/(RTg) ) 


psBcRTs  xi 

'  E  <T,  -  T  )  ^P(-V(CTSI>  <’2> 

L  1  o 


As  a  brief  aside  it  is  interesting  to  note  that,  if  a  constant  char 
temperature  Tg  is  assumed  instead  of  Bqs*  (8),  the  less-accurate  result 
equivalent  to  Eq.  (12)  is 


mGC  =  “  PCBCX1  exp(-Ec/(RTs)) 


(13) 


Equations  (5),  (6)  and  (12)  give 


VP- 


PSBCRTSX1  exP(-Ec/(RTS)) 
Ec(1  -  F  -  r  (PR/Pa,))(T1  -  Tg) 


(14) 


Equation  (14)  may  be  expressed  as 


ln(vp  )  =  In 


p  B  RT  x. 
SCSI 


EC(1  -  p  -  r(PR/P-))(Ti 


Ts)  RTg 


(15) 


The  first  term  on  the  right  hand  side  of  Bq.  (15)  is  likely  to  vary  with  Tg 
much  more  slowly  than  the  second  term,  and  some  experimental  evidence  supports 
this  expectation  [14],  Here,  it  will  be  assumed  that  the  first  term  is 
constant  with  Tg,  at  least  for  a  limited  regime  of  interest.  Thus,  if  the 
variation  of  Tg  with  v  is  measured  (this  would  involve  varying  P, ) ,  a  plot  of 
ln(vpj  against  Tg  should  give  a  straight  line  with  a  slope  equal  to  (-Ec/R). 
This  procedure  is  similar  to  a  thermogravimetric  analysis  115], 
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The  parameter  F  may  then  be  obtained  from  the  intercept  of  the 
straight  line  with  the  axis  Tg^  =  0  if  it  is  assumed  that  ps,  Bq,  x ^  and  T<j 
are  measurable.  This  assumption  is  probably  the  poorest  of  all  those  in  this 
report,  and  it  will  be  further  discussed  here  and  in  Section  4. 


The  three  quantities  ps ,  and  x-j  may  be  obtained  from  post¬ 
ablation  analyses.  This  approach  is  commonly  used  for  ascertaining  heat 
shield  performance  [16].  Although  measurements  of  ps,  the  char  density  at 
the  surface,  and  x^ ,  the  char  layer  thickness,  are  quite  simple  measurements 
for  cold  char,  it  is  necessary  to  assume  that  the  post-ablation  values  are 
indicative  of  the  values  existing  during  ablation.  In  this  respect,  high 
speed  photography  of  the  ablation  at  the  time  that  laser  irradiation  ceases 
may  be  used  to  indicate  if  any  char  loss  continues  after  irradiation. 
Fortunately,  the  thermal  conductivity  of  char  is  usually  higher  than  that  of 
the  uncharred  material  in  zone  B,  and  thus  the  char  layer  thickness  and 
density  may  be  expected  to  "freeze"  during  the  rapid  cooling  that  would  exist 
in  most  environments  of  interest  immediately  following  cessation  of 
irradiation.  The  measurement  of  Bc  involves  thermogravimetric  analysis  of 
the  resulting  char,  and  if  performed  properly  the  analysis  should  yield  a 
value  for  E^,  in  agreement  with  that  obtained  previously.  Such  an  analysis 
should  utilise  temperatures  and  heating  rates  as  high  as  possible,  and  use  an 
inert  atmosphere  in  order  to  simulate  the  convective  blockage  effect  of 
pyrolysis  gases  that  exist  at  the  high  laser-heating  rates.  However,  the 
accuracy  that  may  be  expected  from  this  technique  is  poor  [15]. 


The  quantity  T1  may  be  determined  by  slowly  heating  a  sample  of  the 
resin  in  order  to  find  the  lowest  temperature  at  which  it  is  relatively 
stable.  Typically,  this  temperature  is  only  about  700  K  [13],  and  thus 
T1  <<  Tg.  Hence  errors  in  deducing  are  not  of  major  concern.  As 
mentioned  earlier  in  this  section,  the  behaviour  at  low  heating  rates  may  not 
accurately  represent  that  at  high  heating  rates,  so  some  check  on  the  accuracy 
of  the  above  approach  may  be  warranted.  This  may  be  achieved  by  utilising 
the  frequency  factor  B^  and  the  activation  energy  for  decomposition  of  the 
resin.  If  unavailable,  and  E^  may  be  obtained  by  a  standard 
thermogravimetric  analysis.  Then,  in  analogy  with  the  analysis  leading  to 
Eq.  (12),  mass  continuity  of  resin  for  zone  B  (a  modification  of  Eq.  (4)) 
gives 


vp„  =  +  vp^, 

R  GR  KR1 


=  Mgr  +  (1  -  r)vpR 


pR1BRRT1  ,  ER  . 

~aERr  6XP  (-ir^ 


(16) 
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where  pR  is  the  known  density  of  the  virgin  resin,  pR1  is  the  measurable 
density  of  the  pyrolysed  resin  in  zone  B  near  the  interface  at  x. ,  and 
a  =  dT/dx  in  zone  B  (compare  Bqs.  (8)).  The  quantity  a  is  readily  available 
from  the  heat  conduction  equation  and  the  requirement  for  energy  conservation 
across  the  interface,  i.e. 


V 


(17) 


where  kR  is  the  known  thermal  conductivity  in  zone  B  as  discussed  in  Sec.  2, 
and  kfia  is  equal  to  the  rate  of  heat  flow  into  zone  B;  HR  is  the  enthalpy 
required  to  char  the  resin,  and  is  known  or  easily  measured.  Thus  T1  is 
available  from  Bqs.  (16)  and  (17).  This  completes  the  determination  of  the 
four  quantities  ps,  B^,  xj  and  T-j  shown  by  Eq.  (15)  to  be  necessary  to 
evaluate  F. 


3*4  Determination  of  5 

In  most  cases  with  cw  lasers,  6  will  not  be  small  compared  with 
unity  [ 17,18].  The  simple  gaseous  ablation  products  are  unlikely  to  cause 

absorption,  but  absorption  may  be  significant  in  gases  such  as  C02  and  H20 
formed  from  reaction  of  the  ablation  products  with  the  cold  environmental  gas 
in  front  of  the  ablator  [19].  Experimentally,  these  absorbing  gases  may  be 
precluded  by  using  an  inert  gas  such  as  argon  in  place  of  the  environmental 
gas  of  interest,  thus  making  6=1.  By  measuring  the  ablation  rate  v  in  both 
types  of  atmosphere,  6  for  the  reacting  atmosphere  may  then  be  deduced.  Such 
a  procedure  has  been  performed  in  the  case  of  ablating  samples  of  glass-epoxy 
and  kevlar-epoxy  composites  [7]. 

Alternatively,  6  can  be  obtained  from  Eq.  (2)  and  a  measure  of  v(PL) 
provided  that  PR  is  known  (see  Sec.  3.5).  If  the  results  for  v(PL)  are 
extrapolated  back  to  v  =  0,  then  Eq.  (2)  shows  that  6  =  PR/PL  when  v  =  0. 


3*5  Determination  of  PR 

The  quantity  PR  is  the  remaining  parameter  required  in  Eq*  (2)  in 
order  to  find  Hg.  It  is  limited  by  the  melting  (or  sublimation)  temperatures 
of  ablating  elements.  For  P^  >  30  MW  m~2,  the  surface  temperature  Tg  is 
usually  in  excess  of  3000  K  t5,7],  and  convective  cooling  rates  are  much 
smaller  than  radiative  cooling  rates.  Thus  present  theoretical 
approximations  adequately  determine  the  former  rate  despite  the  large 
uncertainties  [20].  The  radiative  cooling  contribution  is  given 
by  eoT^  where  e  is  the  total  surface  emissivity,  and  a  is  the  Stefan-Boltzman 
constant.  At  4400  K,  oTg  “  20  MW  m“2.  For  ablating  fibre-resin 
composites,  e  >  0.8,  and  thus  if  Tg  is  measured  to  within  200  K,  the  absolute 
uncertainties  in  PR  are  only  around  3  MW  m~2,  which  can  be  ignored  compared 
with  PL  for  PL  >  30  MW  nT2. 
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DISCUSSION  AND  CONCLUSIONS 


At  present  no  analytical  or  numerical  technique,  including  the  one 
described  here,  can  quantitatively  predict  the  significance  of 
thermomechanical  removal  effects  without  some  measurements.  The  analysis 
described  here  shows  how  this  significance  may  be  determined  from  some 
suitable  measurements,  at  least  in  principle,  but  because  some  of  the 
assumptions  in  the  theory  are  poor,  the  main  value  of  the  analysis  is  that  it 
delineates  the  important  physical  and  chemical  parameters  and  clarifies  the 
physical  concepts.  Thus  the  analysis  should  be  considered  to  be  preliminary 
to  more  detailed  work. 


The  model  is  restricted  to  cases  in  which  v  is  observed  to  vary  with 


-2 


Its  main 


P*  in  a  linear  fashion,  and  to  cases  in  which  PL  >  30  MW  m 
shortcomings  are  that  it  assumes  that  the  frequency  factor  Bc  for  the  char  can 
be  determined  experimentally,  and  that  the  char  thickness  can  be  determined 
from  a  post-ablation  measurement.  The  latter  assumption  could  probably  be 
justified  in  some  specific  cases,  but  since  typically  x^  <  2  mm  when 


PL  > 


30  MW  m  *,  the  accuracy  of  post-ablation  measurement  for  x-j  may  even  then 


be  poor.  However,  in  theoretical  treatments  for  the  ablation  of  composite 
materials,  shortcomings  of  this  or  greater  magnitude  are,  by  necessity,  often 
tolerated. 
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FIG.  1  Physical  model  for  the  ablation  of  a  fibre 
composite  material 
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